Condensate fraction and critical temperature of a trapped interacting Bose gas 
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By using a mean field approach, based on the Popov approximation, we calculate the temperature 
dependence of the condensate fraction of an interacting Bose gas confined in an anisotropic harmonic 
trap. For systems interacting with repulsive forces we find a significant decrease of the condensate 
fraction and of the critical temperature with respect to the predictions of the non-interacting model. 
These effects go in the opposite direction compared to the case of a homogeneous gas. An analytic 
result for the shift of the critical temperature holding to first order in the scattering length is also 
derived. 
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The recent experiments on Bose-Einstein condensation 
(BEG) in magnetically trapped atomic vapours ^ have 
stimulated a new interest in the theoretical study of inho- 
mogeneous Bose gases. Although the atom clouds real- 
ized in these experiments are very dilute, the effects due 
^ to the interatomic forces are known to be important at 
low temperature. In particular, the shape and the energy 
' of the condensate cloud ||J^ as well as the dispersion law 
[ — , of the elementary excitations Q are strongly affected by 
' the interaction. In the very recent experiments by the 
Boulder and MIT groups [^-Q , the measured release en- 
ergy and excitation frequencies of the collective modes 
have been found to be in good agreement with the theo- 
retical predictions, thereby revealing important features 
^ , of the trapped Bose condensed gases which are undoubt- 
' edly connected to the interparticle interaction. The ques- 
^1 tion of how two-body forces affect the thermodynamic 
properties of these systems has been also the object of 
several theoretical investigations The critical tem- 
. ^ ' perature of the BEG transition in an homogeneous dilute 
, gas has been recently calculated [|| using the renormal- 
ization group theory. The result is a shift towards higher 
temperatures, with respect to the prediction of the ideal 
Bose gas. Similar results have been also obtained with 
path integral Monte Garlo simulations However, 
no definitive conclusions have been so far drawn on the 
behavior of the condensate fraction, nor of the critical 
temperature in the presence of a confining potential [p] . 
First experimental data on these relevant quantities are 
now becoming available ||5|,|l2[|. 

Finite size effects on the temperature dependence of 
the condensate fraction and on the critical temperature 
in the presence of an external trap have been recently in- 
vestigated by several authors within the non-interacting 
model [fslllj] . In the presence of an anisotropic harmonic 



potential of the form Vext = 'm{u)lx^ -t- w^y^ -I- ujIz'^)/2 
this model predicts, in the large N limit, the well known 
results for the critical temperature 
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where uj = {iUx^^yijJzY^^ is the geometrical average of 
the oscillator frequencies, and for the temperature de- 
pendence of the number of atoms in the condensate 
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Results (0) and (§) are obtained using the semiclassi- 
cal approximation for the excited states and setting the 
chemical potential equal to zero at the transition. The 
first correction to the critical temperature (0) due to the 
finite number of atoms in the trap has been recently 
shown to obey the law 

where uj — [uj^ + lJj, + a;^ ) /3 is the mean frequency. This 
result can be obtained by still employing the semiclas- 
sical description for the excited states, but keeping the 
quantum value /i = i'hu)/2 for the chemical potential at 
the transition. The discretization of the excited energy 
levels gives rise to higher order corrections to 5T^. 

In this Letter we present results for the temperature 
dependence of the condensate fraction and for the crit- 
ical temperature of a dilute Bose gas interacting with 
repulsive forces and confined in an harmonic potential. 
We use a mean field approach and the semiclassical ap- 
proximation for the excited states. This approximation 
is expected to be accurate for temperatures significantly 
larger than the oscillator temperature tiu/kB- This con- 
dition is well satisfied in a useful range of temperatures 
provided N is sufficiently large. In the presence of re- 
pulsive interactions we find that the thermal depletion 
of the condensate is enhanced and the critical tempera- 
ture is shifted towards lower temperatures. These results 
go in the opposite direction compared to the homoge- 
neous case, revealing an interesting behavior exhibited 
by a confined Bose gas. 

Our starting point is the finite temperature generaliza- 
tion of the Gross-Pitaevskii equation within the Popov 
approximation |l5| , [l^ 
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•<?(no(r) + 2nT(r)))<I>(r) = , (4) 



and 



+ 2gn(r) )$'(r,t) + 5no(r)$'*(r, i) . (5) 

The first equation describes the space variations of the 
condensate wave function <I>(r) = {ip{r)) at statistical 
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equilibrium, where ■0(r) is the particle field operator. 
The second one is the equation for the fluctuations of the 
condensate $'(r, i) =< ■ip{r,t) > — <I>(r) which give the 
elementary excitations of the system. InEqs. (^) and (||), 
Vext is the external potential, fj, is the chemical potential, 
g = 47r 7i^a/m is the interaction coupling constant fixed 
by the s-wave scattering length a, hq is the equilibrium 
condensate density no(r) = |$(r)p, n(r) — (V'^(r)-0(r)) 
is the particle density and finally ut is the density of the 
thermally excited particles riT(r) = n{r) — no{r). Equa- 
tions (^), (||) are obtained from the equation of motion for 
the particle field operator ?/;(r,t), by treating the cubic 
interaction term gip'' {r,t)ip{r,t)ip{r,t) in the mean field 
scheme. The Popov approximation consists in neglecting 
the anomalous density mT(i") — (V'(r)V'(r)) — $(rV^ en- 
tering the interaction term. As discussed in Rcf. (l6| this 
approximation is expected to be good for high tempera- 
tures where niT << tit and also in the low-temperature 
regime where tlt and rriT are of the same order, but 
both are negligible. The present mean field approach is 
expected to provide correctly the thermodynamic prop- 
erties of the system, apart from the critical behavior near 
Tc where the mean field approach is known to fail. 

The energies of the elementary excitations can be ex- 
plicitly obtained from Eq. (||) using the semiclassical 
WKB approximation. Let us write the fluctuations of 
the condensate in the form $'(r,t) — Aexp(j(^(r, i)) -|- 

exp(— i(p(r, t)). In WKB approximation the coupled 
equations for $'(r,t) and $'*(r,t), given by Eq. (||) and 
its complex conjugate, can be solved in an explicit way, 
yielding the semiclassical excitation spectrum 

e(P,i-) 

= ^jl^^ + Vext{v)-^i + 2gn{v)^ -^'"^(r), (6) 

where p = KS/ip and e = —hdip/dt are respectively the 
impulse and the energy of the excitation. For a homo- 
geneous system at low temperatures, where nx can be 
neglected, the above excitations coincide with the usual 
Bogoliubov spectrum. The quasi-particles with energies 
e(p, r) are distributed according to the Bose distribution 
function 

exp(e(p,r)/fcBT) - 1 

whereas the particle distribution function can be ob- 
tained from the Bogoliubov canonical transformations 
and is given by 

F{p,r)^-(^) /(p,r) , (8) 
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where the partial derivative is taken at constant conden- 
sate density hq. The non-condensate density is then ob- 
tained by integrating F{p,r) in momentum space 
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and the total number of atoms out of the condensate is 
given by 

Nt = j d^rnrir) . (10) 

In the above equations the chemical potential ^ is fixed 
by the normalization condition 

N^NoiT) + NT , (11) 

where No{T) = J (Pr no(r) is the number of atoms in 
the condensate, with no(r) fixed by the solution of Eq. 

(i- 

We have solved Eqs. (^) - (|l^) in a self-consistent way 
employing the following procedure: 

i) Equation (^) is solved, using the method described in 
Ref. for the condensate density no(r) and the chem- 
ical potential ^, by keeping fixed the number No{T) of 
particles in the condensate and the density ^^(r) of ther- 
mally excited atoms. 

ii) The condensate density and chemical potential found 
in step 1) are used to calculate the excitation energies 
from Eq. (^. A new density nr is then obtained from 
Eq. (^, and new values for the number of atoms out of 
the condensate, Nt, and in the condensate, No{T), are 
derived from Eqs. (10) and (|ll|). 



iii) Step 1) and step 2) are repeated until convergence is 
reached. 

It is worth noticing that the present method accounts 
for finite size effects because of the quantum nature of 
Eq. ^ for the order parameter and of the correspond- 
ing value of the chemical potential. In particular, in the 
absence of interactions, it reproduces result (^ for the 
shift of the critical temperature. 

In Fig. 1 we present our results for the tempera- 
ture dependence of the condensate fraction No{T)/N for 
a system of 5000 Rb atoms interacting with scattering 
length a = llOao, where ao is the Bohr radius, and 
trapped in an axially deformed harmonic potential fixed 
by aHO = (?i/TOw)i/2 = 7.92 x 10"^ cm and by the de- 
formation parameter A — uozl^x — ^z/^y — 's/S. These 
values correspond to the experimental conditions of Ref. 
1^. As clearly emerges from the figure, finite size effects 
are not negligible for this value of TV. However, inter- 
action effects are more important and result in a sizable 
quenching of the condensate fraction. Our results can 
be understood by looking at the role of the two-body 
interaction in Eqs. (^) - (|l^. The chemical potential 
is significantly enhanced, thereby favouring the thermal 
depletion; on the other hand the interaction terms in 
g of Eq. (|) have a minor effect because of the small 
overlap between no(r) and nT(r)- This behaviour differs 
from the one exhibited by a homogeneous dilute Bose 
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gas where one finds a reduction of the thermal deple- 
tion with respect to the prediction of the non-interacting 
model 1,0. 

When the number of atoms in the trap increases, the 
effects due to the interaction become more and more 
important, as explicitly shown in Fig. 2, where we re- 
port results for No(T) corresponding to three values of 
N {N = 10^ 10^ 10^). The calculations of Fig. 2 have 
been carried out using the value 2.56 x lO^'^ for the ratio 
a/cLHO and A = 18/320. This choice corresponds to the 
experimental situation of Ref. for Na atoms. We have 
checked that for such large values of N the results for 
iVo (T) are practically independent of the value of the de- 
formation parameter A. Notice that the quenching effects 
shown in Fig. 2 are quantitatively similar to the ones of 
Fig. 1, despite the much larger values of N contained in 
the trap. This is due to the fact that in the trap of Fig. 
2 (MIT- type trap) the ratio a/ano is a factor 3 smaller 
than in the trap of Fig. 1 (JILA-type trap). 

The shift STc in the critical temperature can be esti- 
mated analytically to the lowest order in the coupling 
constant g by studying the behavior of the trapped gas 
for temperatures T > T^. In this case the order param- 
eter $(r) vanishes and the dispersion of the elementary 
excitations takes the simple Hartree-Fock form 

e(p, r) = + V,^t{r) + 2gnT[v) - ^i . (12) 

To the lowest order in g the total number of particles 
above Tc can be written in the following form 

N = N°{fi) ~2gJ d^r m) , (13) 

where we have used Eqs. (R) - (|l^) with e(p, r) given by 
Eq. (p^). The quantity 72^(r, /z) is the non-condensed 
density given by the non-interacting model 

r4(r, ^i) = A^353/2(exp(-(K.t(r) - ^)/kBT)) , (14) 

where At — fi{2'K /mkBTY^'^ is the thermal wavelength 
and 53/2(2;) = X;,T=i x'^/rv^/'^. 

Bose-Einstein condensation starts at the temperature 
Tc for which the chemical potential reaches the energy 
of the lowest solution of the Schrodinger equation corre- 



sponding to the Haniiltonian (12). To the lowest order 



in g and for large values of N one finds 

/x=^ + 2gnO(r = 0,M = 0) . (15) 

By expanding the r.h.s. of Eq. ( p^ ) around /i = and 
T — T^, we obtain the following relationship for the shift 
6Tc — Tc~ of the critical temperature 



dT 2 J dfi 

-2g / <fv ^ (n° (r = 0) - n° (r)) , (16) 
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where, for simplicity, we have dropped the argument 
/i = from the Uj, functions, here evaluated at T = T°. 
The shift STc of the critical temperature is hence the 
sum of two distinct effects given by the two terms on the 
r.h.s. of Eq. (p^). The first contribution is due to the 
finite number of particles present in the trap and after 
a straightforward integration gives exactly result (||) for 
the shift (5T°/r°. The second contribution instead arises 
from interaction effects and takes the form 

_ 2g J (fr dn^/dn (n^ir = 0) -n^{r)) 
If ~ ~T2 / (Prdn^^/dT 

where S = -,^^-^J^/C(3/2)C(2) ^ 

0.281. By evaluating explicitly the numerical factors one 
gets the relevant result 

/:rpint 

%r - -1-33 ^ N'^' ■ (18) 

The shift 6T^, originating from finite size effects, is al- 
ways negative and vanishes in the large N limit. For 
an axially deformed trap it depends on the deforma- 
tion parameter A = lOz/lOx = ^z/^y through the ratio 
uj/uj (A + 2)/3A^/'^ and is minimum for an isotropic 
trap. Vice versa the shift due to interactions can 

be either negative or positive, depending on the sign of a, 
it increases as N^^^ and does not depend on the deforma- 
tion parameter A, but only on the geometrical average u. 
Furthermore it vanishes for a homogeneous system (see 
Eq. ) where the density rij, is uniform, thereby show- 
ing that a shift in the critical temperature, linear in the 
interaction coupling constant, can occur only in the pres- 
ence of an external trap. It is also worth noticing that in 
the present approach the relationship between the critical 
temperature and the corresponding value of the density 
in the center of the trap is unaffected by the interaction. 
The shift (^8|) is hence always associated with a change 
in the central density produced by interparticle forces. 
The prediction ST^ = (5T° + 6T^"\ obtained from Eqs. 

and (^8|), well agrees with the numerical values ob- 
tained in the self-consistent calculation discussed above 
and reported in Figs. 1 and 2. The relative importance 
of the two effects depends on the value of TV. Notice 
however that the interaction effects depend very weakly 
on N and scale as N^^^a/ano- This behaviour should 
be compared with the effects of the interactions on the 
low-temperature properties of the system (e.g. the size 
and the energy of the condensate) which instead scale ||] 
as Na/ano- 

We finally note that according to the Ginzburg crite- 
rion critical fluctuations violate the mean field result in 
the region 

Sfi^lf^- Mel < ' (19) 
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where /Xj, is the value of the chemical potential at the 
critical point. These fluctuations can affect the shift in 
the chemical potential at the critical point only to order 

. It means that the mean field approach can provide a 
reliable prediction of the shift of the critical temperature 
up to terms linear in the scattering length. 

A more systematic discussion of the thermodynamic 
behavior of a trapped Bose gas, including the tempera- 
ture dependence of the density profiles, the specific heat 
and the superfluid density will be the object of a future 
paper. 

Useful discussions with Franco Dalfovo are acknowl- 
edged. 
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FIG. 1. Temperature dependence of the condensate frac- 
tion for a system of 5000 Rb atoms. The solid line is the 
result for the interacting case, the short-dashed line is the ex- 
act result for the non-interacting case, and the long-dashed 
line corresponds to the non-interacting case in the thermody- 
namic limit (Eq. (2)). 



FIG. 2. Temperature dependence of the condensate frac- 
tion for interacting Na atoms. The solid line corresponds to 
N = 10^, the short-dashed line corresponds to N = 10® and 
the long-dashed line to A'^ = 10^. The dot-dashed line is the 
result for the non-interacting case in the thermodynamic limit 
(Eq. (2)). 
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